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We study within the light-cone path integral approach [3] the effect of the induced gluon 
radiation on high-p^ hadrons in high-energy heavy-ion collisions. The induced gluon spec- 
trum is represented in a new form which is convenient for numerical simulations. For the 
first time, computations are performed with a realistic parametrization of the dipole cross 
section. The results are in reasonable agreement with suppression of high-pr hadrons in 
Au + Au collisions at \fs = 200 GeV observed at RHIC. 

O 

1. One of the most interesting results obtained at RHIC is the suppression of high-pr 
hadrons in Au + Au collisions (for a review of the data, see [1]). It is widely believed that 
parton energy loss due to the induced gluon radiation caused by multiple scattering in 
the quark-gluon plasma (QGP) produced in the initial stage of nucleus-nucleus collisions 
plays a major role in this phenomenon (usually called jet quenching) [2, 3, 4, 5, 6, 7] 
(for a review, see [8]). The most general approach to the induced gluon emission is 
the light-cone path integral (LCPI) approach developed in [3] (see also [9, 10, 8]). It 
accurately treats the mass and finite-size effects, and applies at arbitrary strength of 
the Landau-Pomeranchuk-Migdal (LPM) effect [11, 12]. Other available approaches have 
limited domains of applicability, and can only be used either in the regime of strong (the 
BDMPS formalism [2, 5]) or weak (the GLV formalism [6]) LPM suppression (the GLV 
approach [6], in addition, is restricted to the emission of soft gluons). For this reason they 
can not be used for an accurate analysis of jet quenching for RHIC (and LHC) conditions. 

The LCPI approach expresses the gluon spectrum through the solution of a two- 
dimensional Schrodinger equation with an imaginary potential in the impact parameter 
plane. The imaginary potential is proportional to the cross section of interaction of the 
qqg system (for q — > gq transition) with a particle in the medium, 03 (p) (here p is the 
transverse distance between quark and gluon, the antiquark is located at the center of 
mass of the qg-system). The cr 3 (p) can be written as cr 3 (p) = C(p)p 2 . The factor C(p) 
has a smooth (logarithmic) dependence on p for p <C (hereafter, p D is the Debye 

*Supported by DFG-grant Schi 189/6-1 



1 



screening mass). If one replaces C(p) by C(p e //), where p e ff is the typical value of p, the 
Hamiltonian takes the oscillator form. This approximation, which greatly simplifies the 
calculations, was employed in several analyses [4, 7, 13] (it was also used in the BDMPS 
approach [2, 5]). However, the oscillator approximation turns out to be too crude and 
unsatisfactory. First of all, for RHIC and LHC conditions, the dominating p scale is not 
small enough and the results depend strongly on the choice of p e ff- Another reason why 
the oscillator approximation is unsatisfactory is more serious. In the high energy limit 
the gluon formation length, Lf, becomes larger than the quark pathlength in the QGP, 
and finite-size effects become important. In this regime p e // <C I/^d, and one might 
naively expect that the oscillator approximation should work very well. However, one can 
show [14] that in this regime the dominating N — 1 rescattering contribution (and any 
odd rescattering) evaluated in the oscillator approximation simply vanishes. For RHIC 
and LHC conditions, the finite-size effects play a very important role and the oscillator 
approximation can lead to uncontrolled errors. For this reason, one has to use an accurate 
parametrization of the three-body cross section. It requires numerical calculations for 
solving the Schrodinger equation. 

In the present paper we represent the induced gluon spectrum in a new form which 
is convenient for numerical computations. We, for the first time, calculate the induced 
gluon emission and the nuclear modification factor for RHIC conditions using a realistic 
imaginary potential. 

2. We consider a quark with energy E produced in a medium at z = (we chose the 
z-axis along the quark momentum). The induced gluon spectrum in the gluon fractional 
longitudinal momentum x reads [3] 



dP 

dx 
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2Re J dz l J dz 2 g(x) [K(z 2 , p 2 \z 1 , - K v (z 2 , p 2 \z 1 , p x )\ | =q . (1) 
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Here K is the Green's function for the Hamiltonian (acting in the transverse plane) 
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v(p , z) = .i^Mel, (3) 
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K v (z 2 , p 2 \z x , pj = — — -exp 

2m{z 2 - z\) 



iM(x){p 2 - Pi) 2 i(z2-zi) 



2(z 2 -z 1 ) L f 
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is the Green's function for the Hamiltonian (2) with v{p,z) = 0. In (2), the Schrodinger 
mass is M(x) = Ex(l — x), Lf = 2Ex(l — x)/[m 2 q x' 2 + m 2 g (\ — x)\ is the gluon formation 
length, here m q and m g are the quark and gluon masses that play the role of the infrared 
cutoffs at x ~ 1 and x ~ (in the QGP m q ,g are given by the quark and gluon quasiparticle 
masses). In (3), n(z) is the number density of QGP, and o" 3 is the above mentioned cross 
section of the color singlet qqg system with a particle in the medium. Summation over 
triplet (quark) and octet (gluon) color states is implicit in (3). The a 3 may depend on z 
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(through the Debye screening mass), however, below we will use ^-independent pr>. The 
vertex factor g(x), entering (1), reads 



g(x) 



a s P(x) d d 



2M 2 {x) dp 1 dp 2 ' 



(5) 



where P(x) = Cp[\ + (1 — x) 2 ]/x is the splitting function for the q — > gq transition (Cp 
is the quark Casimir factor). Note that we neglect the spin- flip q — > qg transition, which 
gives a small contribution to the quark energy loss. 

The three-body cross section entering the potential (3) can be written as [15, 3] 



9 1 

^s(p) = gKe(p) + M(! - x )p)\ ~ gM^P) , 



where 



Vqq(p) = o? s C t Cf / dq 



[1 - exp(iqp)] 
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is the dipole cross section for the color singlet qq pair (Ct is the color Casimir for the 
thermal parton (quark or gluon)). 

The spectrum (1) can be rewritten as [L is the quark pathlength in the medium) 



4— = / dzn(z) 
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d^eff /dx (9) can be viewed as an effective Bethe-Heitler cross section, which accounts for 
the LPM and finite-size effects. One can represent (9) as 



d<rfj}(x,z) a s P(x) 
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where the function F is the solution to the radial Schrodinger equation for the azimuthal 
quantum number m = 1 



1 di 



l 



d 



2M(x) \dp / 



. n(z - Qo- 3 (p) 4m 2 - 1 _L 
% 2 8M{x)p 2 L 
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The boundary condition for F(£,p) reads F(£ = 0, p) = s /pa3(p)eKi(ep) , where e = 
[m 2 q x 2 + m 2 (l — x) 2 ] 1 / 2 , and K 1 is the Bessel function. In deriving (10), we used the 
relations [9] 
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The time variable £ in (10), in terms of the variables z and z\ of equation (9), is given by 
£ = z — Z\\ i.e., contrary to the Schrodinger equation for the Green's functions entering 
(1), (10) represents the spectrum through the solution to the Schrodinger equation, which 
describes evolution of the qqg system back in time. It allows one to have a smooth 
boundary condition, which is convenient for numerical calculations. 

3. The jet quenching is usually characterized by the nuclear modification factor (we 
consider the central rapidity region y ~ and suppress the explicit y-dependence) 



da AA {p T )/dydp 2 T 
N bm daPP(p T )/dydp 2 7 



Raa{Pt) = rT jz^TZ wTX2 > ( 12 ) 



where da AA j 1 dydp\ and da pp '/ ' dydp\ are the inclusive cross section for A + A and p + p 
collisions, and N bin is the number of the binary nucleon-nucleon collisions. The effect of the 
parton energy loss on the high-p^ hadron production in A+A collisions can approximately 
be described in terms of effective hard partonic cross sections, which account for the 
induced gluon emission [16]. Using the power-low parametrization for cross section of 
quark production in p + p collisions oc l/pp PT ^ one can obtain 

1 1 7 J / \ 

R AA (p T ) « P (p T ) + ——J dzz n ^- 2 D h q (z,^) J dx(l-x) n <rr/*)-* ^ , (13) 
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J(p T )= Jdzz n ^- 2 D h q (z,^), (14) 



where P is the probability of quark propagation without induced gluon emission, 
dI(x,pT)/dx is the probability distribution in the quark energy loss for a quark with 
E = pj<, D q (z,pT/z) is the quark fragmentation function. Note that, since n(pT) ^> 1, 
the z-integrands in (13), (14) are sharply peaked at z ~ z (z is the value of z at which 
z n ( pT ~ 2 Dq(z, Pt/ z) has a maximum). For this reason (13) to quite good accuracy can be 
approximated as 



1 dlix PT ) 

Raa( P t) « Po(Pt) + J dx(l - x)<ptI-*)- 2 . 



We take the Pq and spectrum in the radiated energy entering (13) in the form 
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where x m i n ~ rn g /E, and it is assumed that the spectrum equals zero at x < x m i n - For- 
mula (17) corresponds to the leading order term of the series in L/L ra( i (here L ra d is the 
radiation length corresponding to the energy loss AE ~ E) of the spectrum derived in [17] 
for the photon emission and, strictly speaking, is only valid for AE <C E. However, even 
for the more broad domain AE < E (which is interesting from the point of view of jet 
quenching at RHIC) (17) reproduces the energy loss spectrum evaluated assuming inde- 
pendent gluon radiation to a reasonable accuracy. An accurate calculation of the nuclear 
modification factor accounting for the higher order terms in L/L ra d in the approximation 
of independent gluon emission [16] does not make sense because this approximation it- 
self does not have any theoretical justification. Note that our spectrum is automatically 
normalized to unity. 

The effective exponent u{pt) for quark production entering (13) is close to that for 
hadron production, nh{pr)- A small difference between these quantities (stemming from 
the pt dependence of the integral (14)) is given by u(pt) — n^Pr) = din J(p T )/dliap T . 
For Au + Au collisions at yfs = 200 GeV we use rihipr) = npr/ipr + b) with n = 9.99 
and b = 1.219 corresponding to the parametrization dajdydp\ = A/(p T + b) n obtained 
in [18] for 7r° production in p + p collisions. The above procedure allows one to avoid 
uncertainties of the pQCD calculations of the partonic cross sections. 

4. To fix the fn qt g and fi£> we use the results of the analysis of the lattice calculations 
within the quasiparticle model [19] . For the relevant range of temperature of the plasma 
T ~ (1 — 3)T C (T c pa 170 MeV is the temperature of the deconfinement phase transition) 
the analysis [19] gives for the quark and gluon quasiparticle masses m q pa 0.3 and m g pa 
0.4 GeV. With the above value of m g from the perturbative relation ji D — y/2m g one 
obtains fio ~ 0.57 GeV. To study the infrared sensitivity of our results we also perform 
computations for m g = 0.75 GeV (with /i£> = 0.57 GeV). This value of the infrared cutoff 
for gluon emission in parton-nucleon interaction has been obtained from the analysis of 
the low-rr proton structure function within the dipole BFKL equation [15, 20]. It seems 
to be reasonable for gluon emission in the developed mixed phase and for fast gluons with 
Lf > L. This value agrees well with the natural infrared cutoff for gluon emission in 
the vacuum m g ~ l/R c , where R c pa 0.27 fm is the gluon correlation radius in the QCD 
vacuum [21]. The above two values of m g give reasonable lower and upper limits of the 
infrared cutoff for the induced gluon emission for RHIC and LHC conditions. 

We perform numerical calculations for fixed and running a s . In the first case we take 
a s = 0.5 for gluon emission from light quarks and gluons, and a s = 0.4 for the case of 
c-quark. For the running a s , we use the parametrization (with Aqcd = 0.3 GeV) frozen 
at a s = 0.7 at low momenta. This parametrization is consistent with the integral of a s (Q) 
in the interval < Q < 2 GeV obtained from the analysis of the heavy quark energy loss 
(in vacuum) [22]. Previously this parametrization was used to describe successfully the 
HERA data on the low-x proton structure function within the dipole BFKL approach 
[15, 20]. To incorporate the running a s in our formalism, we include a s in the integrand 
on the right-hand side of (10) and take for virtuality Q 2 = aM{x)/$ > . The parameter a 
was adjusted to reproduce the N — 1 rescattering contribution evaluated in the ordinary 
diagrammatic approach [23]. It gives a pa 1.85. For the dipole cross section (7) we take 
Q 2 = q 2 . 
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We assume the Bjorken [24] longitudinal expansion of the QGP with T 3 r = T 3 r which 
gives n(z) oc \jz for z > r . We use the initial conditions suggested in [25]: T = 446 
MeV and r = 0.147 fm for RHIC, and T = 897 MeV and r = 0.073 fm for LHC. 
For RHIC, the above condition were obtained from the charged particle pseudorapidity 
density dN/dy 1260 measured by the PHOBOS experiment [26] in Au + Au collisions at 
y/s = 200 GeV assuming an isentropic expansion and rapid thermolization at r ~ 1/3T . 
The LHC parameters correspond to dN/dy ~ 5625 at = 5.5 TeV, which was estimated 
in [27]. The above initial conditions for RHIC (translated into r = 0.6 fm) agree with 
those used in successful hydrodynamic description of Au + Au collisions at RHIC [28]. 
Note that, since the dominating p-scale in (9) oc \fz for z <^ Lf, our results are not very 
sensitive to r (for a given entropy). The maximum parton pathlength in the hot QCD 
medium is restricted by the life-time of the QGP (and mixed) phase, r max . 1 . We take 
Tmax ~ Ra ~ 6 fm. This seems to be a reasonable value for central heavy-ion collisions, 
since, due to the transverse expansion, the hot QCD matter should cool quickly at r > Ra 
[24]. 

In Fig. 1 we show the induced gluon spectrum for q — > gq transition for RHIC condi- 
tions for the quark pathlength L = 6 fm obtained with m q = 0.3 and m q = 1.5 GeV for 
m g = 0.4 and m g = 0.75 GeV. In Fig. 1 we also show the Bethe-Heitler spectrum (dashed 
line). One sees that the LPM and finite-size effects strongly suppress the gluon emission. 
The gluon emission from c-quark is suppressed in comparison with light quark due to 
larger mass, which leads to decreasing of the dominating p scale (note that the spectrum 
is not sensitive to the light quark mass, except for x ~ 1). One can see from Fig. 1 that, 
although the Bethe-Heitler spectrum differs strongly for two values of m g , the difference 
becomes relatively small for the spectrum, which accounts for the LPM and finite-size 
effects. It is connected with the fact that, due to the multiple scattering and finite-size 
effects, the dominating 1/p-scale becomes larger than m g ; namely this in-medium scale 
plays the role of the infrared cutoff at high energies [4] (however, of course, for not very 
high pt the value of m g is still important). We do not show the spectra for running a s . 
They are similar in form (but somewhat suppressed at moderate fractional momenta). 
The LPM suppression for LHC is considerably stronger than for RHIC, but the spectra 
are similar in form, and we do not show them as well. 

In Fig. 2 we plot the quark energy loss AE = E f£ . dxxdP/dx evaluated for fixed 
(solid line) and running (dashed line) a s for RHIC and LHC conditions for L = 6 fm with 
m g = 0.4 GeV (thick lines) and m g = 0.75 GeV (thin lines). The results for a s = 0.5 
agree roughly with that for running a s for E < 10 GeV but at higher energies the energy 
dependence is steeper for fixed a s . This says that the typical p-scale becomes smaller 
with increasing energy. It is also seen from the decrease of the relative difference between 
the curves for m g = 0.4 and 0.75 GeV. 

In Fig. 3a we compare the nuclear modification factor (13) for To = 446 MeV calculated 
using the NLO KKP fragmentation functions [29] for running a s with that obtained 
at RHIC [1] for central Au + Au collisions at y/s = 200 GeV. For illustration of the 

1 For our choice of the initial conditions the life-time of QGP is <~ 3 fm for RHIC. However, in the 
interval r ~ 3 — 6 fm the density of the mixed phase is practically the same as that for the pure QGP 
phase. 



6 



dependence on T , in Fig. 3b we also present the results for T = 375 MeV. The theoretical 
curves were obtained for L = 4.9 fm. It is the typical parton pathlength in the QGP (and 
mixed) phase for r max = 6 fm. We present the results for m g = 0.4 and 0.75 GeV. For 
Pt ^ 15 the results for a s = 0.5 are close to that for running a s and we do not plot them. 
The results for the quark (solid line) and gluon (dashed line) jets are shown separately 
(note that for y/s = 200 GeV the quark and gluon contributions are comparable). The 
suppression is somewhat stronger for gluon jets. One can see from Fig. 3a that the 
theoretical Raa for m g = 0.4 is in reasonable agreement with the experimental one. One 
should bear in mind, however, that our calculations neglect the collisional energy loss 
[30]. For pt ~ 5 — 15 GeV the collisional energy loss may increase the total energy loss by 
~ 30 — 40%. In this case (if one takes the initial conditions [25]) the value m g = 0.75 GeV 
would be more preferable for agreement with the RHIC data. As mentioned previously, 
this value is reasonable for the mixed phase and for gluons with Lf > L. Since, for 
y/s = 200 GeV, the medium spends about half of its time in the mixed phase, the 
effective infrared cutoff may be larger than the gluon quasiparticle mass in the QGP. 
For this reason, the collisional energy loss may be included without using an unrealistic 
infrared cutoff for the induced energy loss. The possible remaining small disagreement 
with the data may be avoided by taking a somewhat smaller value of T (or a s ). In any 
case, it is clear that, for such a complicated phenomenon, it is hardly possible to expect 
a perfect agreement with experiment and the agreement found in the present paper is 
surprisingly good. 

The above estimate for the collisional energy loss has been obtained for the pQCD 
plasma. Presently, there are some indications [31] that the medium produced at RHIC 
may be a strongly coupled QGP. The radiative energy loss should not be very sensitive to 
the dynamics of the QGP (for the same number density of the QGP). However, it may be 
important for the collisional energy loss. Unfortunately, the corresponding calculations 
have not been made yet. It is interesting that our results give support for the scenario with 
strongly coupled QGP. Indeed, this scenario requires a s > 0.5 [31] for thermal partons. 
The Raa is sensitive to the radiation of soft gluons with an energy of about a few \i D . 
One can expect that, for such gluons, a s should be close to that for thermal partons. 
We obtained agreement with the data with a s , which is frozen at a value of 0.7 at low 
momenta. If a s is frozen at a value below 0.4-0.5, the theoretical Raa strongly disagrees 
with that observed at RHIC. 

5. In summary, we have represented, within the LCPI approach [3], the induced gluon 
spectrum in a new form convenient for numerical calculations and carried out computa- 
tions of the induced gluon emission from fast partons in the expanding QGP for RHIC 
and LHC conditions. The calculations for, the first time, have been performed with a 
realistic parametrization of the dipole cross section. The theoretical nuclear modification 
factor calculated for the initial conditions obtained from the charged particle rapidity 
density observed at RHIC [25] and the hydrodynamic simulation of the RHIC data [28] 
is in a reasonable agreement with that observed at RHIC. 
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Figure 1: The induced gluon spectrum (solid line) for q — > gq transition versus the gluon 
momentum k for RHIC conditions for E — 5, 10 and 20 GeV, L = 6 fm obtained using 
(8), (10) with a s = 0.5 for m q = 0.3 GeV (a,b) and m q = 1.5 GeV (c,d); m g = 0.4 GeV 
(left) and 0.75 GeV (right). The Bethe-Heitler spectrum is shown by the dashed line. 
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Figure 2: The energy dependence of the quark energy loss for RHIC (a) and LHC (b) for 
L = 6 fm obtained with a s = 0.5 (solid line) and running a s (dashed line), m g = 0.4 GeV 
(thick lines) and m g = 0.75 GeV (thin lines), m q = 0.3 GeV. 
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Figure 3: The nuclear modification factor (13) for central Au + Au collisions at ^/s = 200 
GeV for quark (solid line) and gluon (dashed line) jets obtained with m g = OA GeV 
(thick lines) and m g = 0.75 GeV (thin lines) for running a s . The experimental points 
(from [1]) are for the following: circle - Au + Au -> tt° + X (0-10% central) [PHENIX 
Collaboration], square - Au + Au -> + X (0-10% central) [PHENIX Collaboration], 
star - Au + Au ^ h ± + X (0-5% central) [STAR Collaboration]. 



12 



